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5 ■ Abstract 

o 

Using ab-initio computational techniques on crystal determined 
clusters, we report on the similarities and differences of Al5o(C5(CH3)5)i2, 
Ga23(N(Si(CH3)3)2)n, and Auic^SCyC^Hs^ ligand-protected enis- 
le ■ ters. Each of the ligand-protected clusters in this study show the 

*n ' similar stable character which can be described via a electronic shell 

00 

model. We show here that the same type of analysis leads consistently 
to derive a superatomic electronic counting rule, independently of the 
<^ , metal and ligand compositions. One can define the cluster core as the 

set of atoms where delocalized single-angular-momentum-character or- 
bitals have hight weight using a combination of Bader analysis and 
the evaluation of Khon-Sham orbitals. Subsequently one can derive 
^ | the nature of the ligand-core interaction. These results yield further 

insight into the superatom analogy for the class of ligand-protected 
metal clusters. 
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1 Introduction 



Protected gold clusters have been synthesized in several sizes and compo- 
sitions. Due to their intrinsic stability and potential applications in nan- 
otechnology, they have received broad interest in past decades. For example, 
size dependent optical, electrochemical and catalytical properties have been 
experimentally determined [IJISIEIIIIEIEIITIIE]. In addition, these clus- 
ters can be expected to be organized as building blocks of materials with 
new interesting properties. Recent breakthroughs such as determining the 
crystal structure of ligand-protected clusters [9] [101 HH [121 US] , has led to fur- 
ther investigations of their electronic structure via ab-initio simulations . As 
a result, the electronic structure and derived properties of ligand-protected 
gold clusters can be modeled via a modified electronic shell model, termed 
superatom model [HI [151 EE E] • 

Electronic shell models have been successfully used to understand and 
predict properties of bare metallic clustersJTSJIISlEn]- One possible spherical 
shell model is a 3D harmonic oscillator with an anharmonic, angular momen- 
tum dependent, term. From this model several properties are derived, most 
notably the high stability of some clusters. For a given composition, the 
model predicts stable clusters with large HOMO-LUMO gaps correspond- 
ing to the closing of an electronic shell level. The order of the shells and 
the magic numbers in this spherical shell model depends then on the anhar- 
monic parameter. For an anharmonic parameter in the intermediate region 
the order of the shells is: 

IS 2 IP 6 ID 10 2S 2 IF 14 2P 6 1G 18 2D 10 3S 2 1H 22 2F U 3P 6 1J 26 2G 18 

predicting magic numbers at 58 and 138 electrons, among others, correspond- 
ing to the closing of the 1G and II shell respectively (U = 0.03, equation 

Aluminum and gallium metalloid clusters (clusters containing more metal- 
metal bonds than metal-ligand bonds) have also been characterized experi- 
mentally and theoretically in an effort to understand how properties evolve 
from clusters to bulk and their stability has been explained through various 
models [2T| [221 [23]. However, reports on bare aluminum clusters have shown 
the most stable species have superatomic character with a magic number of 
electrons, which adheres to the electronic shell model [2~4"1 125]. For example, 
the AlJ~ 3 cluster is resistive to O2 etching with 40 electrons (magic number), 
while its neutral counterpart (Al 13 ) is defined as a superhalogen based on its 



2 



similar electron affinity to halogens on the periodic table [2U ETJ ESI EH]- 
This result, along with the studies on ligand-protected gold clusters, suggest 
the electronic shell model could describe metalloid clusters composed of alu- 
minum and gallium. Recently our group has analyzed aluminum metalloids 
and successfully illustrated a superatom electronic structure exists which re- 
lates to a cluster's overall stability [SU]. For few gallium metalloid clusters, 
it has been predicted the electronic shell model may be successful; however 
to our knowledge no theoretical investigations into the electronic structure 
have been performed. 

In order to derive an electronic shell structure for nearly spherical ligand- 
protected clusters, we currently use an orbital projection on spherical har- 
monics integrated on the cluster region[14j. The orbitals can be labelled with 
a given angular momentum using the coefficients and one can determine the 
order of electronic shells in the cluster system. It is necessary however to 
take the role of the ligand layer into account to determine the expected 
number of delocalized orbitals participating in the shell structure. In the 
superatom model, one takes the ligand-core interaction into account. In a 
ligand-protected superatom cluster the core atoms participate collectively to 
give rise to a cluster size delocalized orbitals; while the protective units par- 
ticipate only in localized or interface bonding states. The protective ligands 
can behave by either depleting or donating electrons (or remaining neutral) to 
the superatom electronic structure and subsequently will give rise to different 
electron counts i.e-magic numbers. The core-ligand interaction is not only 
important from a theoretical perspective, but proves to have implications in 
experimental observations as well. [3TJ [17] 

Completing the characterization of the clusters as ligand-protected super- 
atoms requires an understanding of not only the electronic structure, but the 
ligand-core interaction as well. We show using ab-initio computational tech- 
niques, three different ligand-protected clusters (Al 50 Cp^ 2 Cp*=C5(Cn 3 ) 5 
and Ga23Ln, L = N(Si(CH 3 )3) 2 and Auio2(SR)44, R=C702H 5 ) can be fully 
characterized as superatoms. Projection of the Kohn-Sham (KS) orbitals on 
a local atomic basis is used to determine which atomic layers contribute to 
form the delocalized orbitals and are part of the cluster core. Using Bader 
analysis we find one can characterize the local atomic electronic structure and 
gain information into the nature and description of the ligand-protecting shell 
for each cluster. Based on these analyses of the Mn[Xy] z clusters, one can 
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predict the number of delocalized electrons using the equation: 



n * = NV M -YV x -z 



(1) 



where N is the number of atoms in the core with valence Vm, Y is the 
number of protective units depleting V x electrons each and the cluster has 
an overall charge of z, that contribute to the superatomic orbitals in the 
superatom model. Finally, we compare the similarities and differences of the 
three superatom complexes and their expected properties. 

2 Computational methods 

The computations were done using the GPAW code, which performs calcula- 
tions based on Density Functional Theory [321 The code is a grid based 
implementation of the projector- augmented wave method (PAW). Further- 
more a frozen core approximation is used. H(ls), C(2s2p), Al(3s3p), S(3s3p), 
0(2s2p), Au(5dls), N(2s,2p), Si(3s,3p), Ga(4s,4p) are treated in the valence. 
The exchange-correlation functional used in all calculations is PBE [31] . Re- 
laxation of the system is performed until the forces around all atoms are 
below 0.05 eV/A. The crystal structures for the Aui 2(SR)44, Al 50 CpJ 2 and 
Ga23Ln were obtained from the experimentally reported structures through 
the CCDC database[35j [361 [9J- From the crystal structure the coordinates 
of a single cluster were isolated and were allowed "in vacuum" to optimize 
without constraints. 

Atomic charge state (extra or missing local charge with respect to the 
atomic number) is determined using a Bader type of analysis [37J EH] • Pro- 
jecting the all electron partial waves (i.e. the wave functions of the isolated 
atoms) into molecular orbitals (PLDOS) is done within the PAW formalism 
following [33] . 

The superatomic analysis is done by a projection of the KS orbitals on 
spherical harmonics [2]. For simplicity the origin is chosen in the center 
of mass of the cluster. For a given core radius Ro, the angular momentum 
weight q associated to the KS orbital ip is defined using: 





m=—l 



rRo r 



fr= / Yr*(e,<pMr,9,<p)dn r 2 dr. 



(3) 



Jo Jo. 
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where Y™{0, (p) is a spherical harmonic function with degree / and order m 
and dQ = sin 9d(pd9. 

3 Results and discussion 

3.1 Protected gold clusters Auio2(SR)44 

First it is important to see in the protected gold cluster that the gold atoms 
can be divided in two sets: the core and the ligand set [39]. As this result 
has been derived previously and generalized to other protected gold clus- 
ters we include here the discussion for completeness and as an ilustration 
of the method that will be utilized in the aluminum and gallium case. The 
Auio2(SR)44 cluster can be viewed in subsequent layers, Figured] (a). We 
find from the analysis that the atoms in the first three layers r\ to r 3 have 
each a very small average Bader charge, 0.01, 0.00 and 0.06 |e| respectively in 
Figure [2] (a). The atoms in the outermost layer show a small but distinct 
positive 0.13 |e| mean charge. The projection on an atomic basis for the 
inner gold atoms (layers r\ to r 3 ) shows a s-p hybrid band with high weights 
around higest occupied molecular orbital (HOMO) in the PLDOS Figure [3] 
(a). In contrast there are almost no weight on the orbitals around HOMO 
in the external layer r&. Those KS orbitals also have high s-p atomic local 
component, are delocalized superatom orbitals. The superatomic projection 
shows a change of angular momentum G to H and corresponds to a clear 
gap of 0.48 between the HOMO and the lowest unoccupied molecular orbital 
(LUMO) [H], see Figure (a). As a consequence, the superatom core is 
composed of the first three shells r\ to r 3 only; the core being covered by 21 
units (Au(SR)2)i9 and (Au2(SR) 3 )2 including the gold atoms belonging to 
the shell r 4 . One can rewrite then the cluster formula, making explicit the 
core-ligand structure as follows: Au7 9 [(Au(SR)2)i9(Au 2 (SR) 3 )2] implying 79 
gold atoms contribute with one electron (6s) to the delocalized super atom 
counting but 21 of those are depleted by the protective units. This gives 
there are 58 superatom electrons, corresponding to a cluster with electronic 
closed shell configuration in a spherical potential. 
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3.2 Protected aluminum Al 5 o(Cp*)i2 

The aluminum metalloid cluster can be view geometrically into three distinct 
shells of aluminum atoms, an inner Alg layer, encapsulated by 30 Al atoms 
with an exterior shell of twelve Al atoms Figured] (b). The exterior twelve 
aluminum atoms are bonded to twelve pentamethylcyclopentadienyl (Cp*) 
ligands. It should be noted, whether the exterior layer should be considered 
part of the ligand or the core is currently a subject of discussion J23JH0]. In 
order to understand if the outer twelve aluminum atoms are part of the ligand 
or not, we performed a similar ab-initio simulation analysis to determine the 
core and ligand sets as well as the electronic state of the atoms. 

The twelve exterior Al atoms have an average Bader charge of 0.91 |e| 
(Figure[2] (b)), while the innermost atoms has zero or negative charge. Thus, 
these 12 Al atoms donate 1 electron to the electron withdrawing ligand Cp*. 
The projection of KS orbitals on local atomic basis is described in Figure [3] 
(b) where it is shown that all shells, r x to r 3 , contribute to the delocalized 
superatomic orbitals. The outermost shell near the HOMO level shows 
more s states than p per atom, which can be attributed to the loss of 1 p- 
electron per Al atom and corresponds to the Bader charge results. Further, 
the projection on spherical harmonics shows the Al 50 (Cp*)i2 has an electronic 
structure corresponding to the one expected for a superatom cluster with a 
shell closing of 138 electrons as reported previously [SU] see Figure [5^b). The 
HOMO-LUMO gap (0.94 eV) corresponds to the transition from the II shell 
to the 2G shell. The splitting of those II and 2G shells can be obtained 
within the spherical shell model with the anharmonic parameter indicated in 
the introduction. Further, when changing the radius parameter in the pro- 
jection of spherical harmonics we confirm it is important to include the most 
outer aluminum atoms since they contribute to form the higher delocalized 
superatomic orbitals (Figure 1 in Supporting Information). Thus, the combi- 
nation of these results point to the characterization of this metalloid cluster 
as a ligand-protected superatom with ionic bonding to the Cp* ligands. The 
structural separation in core-ligand can be represented writing the chemical 
formula as Al 50 [Cp*i2]. The 50 Al atoms of the core contribute each with 
three electrons while each protective ligand Cp* deplete one electron giving 
a total number of superatom electrons of 138. 



6 



3.3 Protected gallium Ga23Rn 

The Ga23(N(SiMe3)2)n metalloid cluster can also be viewed as a cluster 
containing multiple radial atomic layers as in the previous two cases [1] (c). 
The inner most layer consists of one gallium atom followed by two consecutive 
layers of eleven gallium atoms, with the outer eleven being surrounded by 
ligands. Previously, the gallium metalloid cluster has been described as a 
Gai2 core surrounded by 11 GaL units [35]. We find the 11 Ga atoms on the 
exterior layer have an average Bader charge of 0.52 |e|, with the inner shell 
having a negative -0.40 and neutral 0.04 |e| mean Bader charge (Figure 12(c)). 
The Bader charge value is indicative of polarized bonding, which should come 
as no surprise since the ligand is composed of nitrogen. It is well-known that 
the nitrogen atom participate in polarized bonding. From the projection on 
local atomic basis PLDOS (Figure |3] (c) ) it is obtained all gallium atoms 
are in the same electronic state and contribute to the delocalized superatom 
states around HOMO. The combination of Bader analysis and PLDOS allows 
us to conclude the gallium cluster should be described as a metallic core 
of 23 atoms protected by 11 electron depleting units. The projection of 
the KS orbitals on spherical harmonics reveals the ligand-protected gallium 
cluster does adhere to the superatom model with a closed 1G shell with 
58 superatomic electrons (Figure [5] (c)). The cluster has a large gap of 
1.34 eV, which is indicative of its stable nature as well. The derived core- 
ligand composition is thus Ga23[Rn] with 23 gallium atom contributing with 
3 electrons each and 11 protective units depleting one electron and thus 
corresponding to a counting of 58 superatom electrons. 

3.4 Comparison 

All three clusters considered in this study have a core that can be separated 
in concentric layers. Each layer being formed by metal atoms having not 
only the same radius but also the same local Bader charge. Likewise, their 
electronic structure can be explained using the superatom model with the 
number of electrons contributing to the delocalized superatomic orbitals ob- 
tained using equation [TJ Using this electron counting rule one obtains 58, 
138 and 58 electrons for the Auio2(SR)44, Al 50 [Cp*i2 and Ga23(N(SiMe3)2)n 
clusters, respectively. Each of the cluster are superatom complexes, but the 
projection on the local atomic basis shows differences in the layers that con- 
tribute to the delocalized superatomic orbitals. For the gold cluster, the gold 
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M N [X Y ] 



n e = NV M -YV x -z 



Au 79 [(Au(SR) 2 ) 19 (Au 2 (SR) 3 ) 2 ] 

Al 50 [Cp* 12 ] 

Ga 23 [Rn] 



n e = 79 - 21 = 58 

n e = 3 x 50 - 12 = 138 

n e = 3 x 23 - 11 = 58 



Table 1: Formulas and derived superatomic counting rule 



atoms in the protective exterior do not participate to form such orbitals. In 
the cases of the metalloid superatom complexes their exterior metal atoms 
do contribute to form the superatomic orbitals. Further differences can be 
seen from the Bader analysis. The charge distribution in the ligand-protected 
aluminum and gallium clusters proceeds from negative to neutral to positive 
from the center to outer shells. However, in the ligand-protected gold cluster 
most inner two shells (ri and r 2 ) are neutral, with the third shell r 3 being only 
slightly positive and the exterior shell r 4 positive as well. Another difference 
is found when focusing on the ligand and exterior atomic layer interaction. 
The exterior metal atoms in all three cases have very different Bader values 
with the aluminum value being the largest (0.91 |e|) and gold the smallest. 
It is interesting that the gallium value is 0.52 |e|, yet still positive as for both 
aluminum and gold. These differences point to the nature of the ligand-metal 
bonding in the superatom complexes. It is well-known that the bonding for 
many ligand-protected gold clusters is covalent. We have shown one can have 
ionic and polarized covalent bonding in other superatom complexes as seen 
in the aluminum and gallium cases presented here. 

4 Conclusion 

We have investigated three distinct ligand-protected clusters composed of 
gold, aluminum and gallium. The electronic structure of these clusters can 
be described via the superatom model and have magic numbers correspond- 
ing to 58, 138 and 58 electrons, respectively. One can also define the metallic 
core as the group of atoms which participate in the delocalized superatomic 
orbitals using the analyses presented here. Finally, we present evidence that 
the charge distribution in the core varies strongly with the clusters and do 
not affect the validity of the electronic shell model. The charge variation is 
related to the metal-ligand bonding that also varies with the clusters. We 
found a covalent, ionic and polarized covalent bond for the Auiq 2 (SR)44, 
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Al 5 oCp*i2 and Ga 2 3L n clusters respectively. As previously and succesfully 
done for the ligand-protected gold clusters, consequences into optical and 
charging properties, voltammetry, NMR experiments and reactivity should 
be expected for the both the aluminum and gallium (as well as other com- 
postions) superatom complexes. We hope this study will lead to further 
theoretical and experimental investigations. 
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Figure 1: Radial atomic decomposition of the ligand protected clusters con- 
sidered. Metal atoms are represented with spheres and all others with sticks. 
Metal atoms in the same radial layer share the same color. Some atoms 
have been removed to enhance the view of the core, (a) Protected gold clus- 
ter Auio2(SC 7 2 H5)44: ri, r 2 , r 3 and r 4 are colored in yellow, pink, blue 
and violet respectively, (b) Protected aluminum cluster A1 50 (C 5 (CH3)5) 1 2: 
r\ and r 2 are colored in violet and r 3 in blue, (c) Protected gallium clus- 
ter Ga23(N(Si(CH3) 3 ) 2 )n: ri, r 2 , r 3 are colored in yellow blue and violet 
respectively. 
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Figure 2: Bader charge as a function of the radial position of the atoms. 
Negative charge indicates excess with respect to the neutral atom charge. 
The origin of the coordinates is chosen in the center of mass of the cluster. 
Metal atoms are all included individually. The rest of the charge is indicated 
as a per ligand average, (a) Protected gold cluster Auit^SCyC^Hs^ (b) 
Protected aluminum cluster Al5 (C5(CH3) 5 )i 2 (c) Protected gallium cluster 
Ga 23 (N(Si(CH 3 )3)2)ii. 
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Figure 3: Projected Local Density of States using atomic basis. The pro- 
jection is integrated in the radial layers composed by metal atoms, (a) 
Protected gold cluster Auh^SC/C^Hs^ (b) Protected aluminum cluster 
Al 5 o(C5(CH 3 )5)i2 (c) Protected gallium cluster Ga 2 3(N(Si(CH 3 ) 3 ) 2 )ii. 
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Figure 4: Super-atom analysis: the angular momentum coefficient q from 
Eq [2] (/ = 0,1. ..6) as a function of the energy of the projected KS or- 
bital. The difference between the sum of the coefficients and the norm of 
the KS orbital inside the sphere of radius Ro is A. (a) Protected gold cluster 
Auio2(SC70 2 H 5 )44, the value of i? is 7.5 A(b) Protected aluminum clus- 
ter Al 5 o(C5(CH3)5) 12 , the value of R is 6 A(c) Protected gallium cluster 
Ga 2 3(N(Si(CH 3 )3) 2 )i 1 , the value of R is 5.5 A. 
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Figure 5: Supporting Information figure. Super-atom analysis: the angular 
momentum coefficient q from Eq 1(7 = 0, 1...6) as a function of the energy of 
the projected KS orbital. The difference between the sum of the coefficients 
and the norm of the KS orbital inside the sphere of radius Rq is A. The 
analysis is shown for the protected aluminum cluster A1 50 (C5(CH 3 )5) 1 2 with 
R core radius fixed at (a) 3 A, (b) 5 Aand (c) 6 Awich corresponds to the 
half distance between radial layers. 
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